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ABSTRACT 

Higgs singlet superfields are usually present in most extensions of the Minimal Supersym- 
metric Standard Model (MSSM) that address the /i-problem, such as the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM) and the Minimal Nonminimal Supersymmetric 
Standard Model (MNSSM). Employing a gauge- and flavour- covariant effective Lagrangian 
formalism, we show how the singlet Higgs bosons of such theories can have significant con- 
tributions to 5-meson flavour- changing neutral current (FCNC) observables for large values 
of tan (3 > 50 at the 1-loop level. Illustrative results are presented including effects on the 
Bg and B^ mass differences and on the rare decay Bg fi~^fi~ . In particular, we find that 
depending on the actual value of the lightest singlet pseudoscalar mass in the NMSSM, the 
branching ratio for Bs —>■ can be enhanced or even suppressed with respect to the 

Standard Model prediction by more than one order of magnitude. 
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1 Introduction 



Supersymmetry (SUSY), softly broken at the electroweak (EW) or TeV scale, provides one 
of the most self-consistent frameworks that enables one to technically address the so-called 
gauge hierarchy problem. In particular, SUSY protects the EW scale from quantum correc- 
tions that could be induced by possible new dynamics at the Grand Unified Theory (GUT) 
or Planck scale. Whilst the exact mechanism of SUSY breaking remains unknown so far, 
supersymmetric theories do in general contain many sources of FCNC's amongst the soft 
SUSY-breaking terms. In order to maintain agreement with experimental data, the hy- 
pothesis of Minimal Flavour Violation (MFV) is frequently assumed, where the well-known 
Glashow-Iliopoulos-Maiani (GIM) mechanism can be naturally implemented [1]. Within 
the MFV framework, all phenomena of flavour and CP violation in the theory originate 
from the Cabbibo-Kobayashi-Maskawa (CKM) mixing matrix V [2], such that all FCNC 
interactions vanish in the limit V — > I3. 

The FCNC effects on B- and K-meson observables mediated by Higgs bosons have 
been well studied in the MSSM [3-7]. At large values of tan/3 = V2/V1, the ratio of the 
Higgs doublet vacuum expectation values (VEVs), the couplings of the Higgs bosons to 
down-type quarks becomes enhanced and 1-loop threshold effects can make a significant 
contribution to FCNC processes, such as the Bg — Bg and B^ — Bd mixings, and the leptonic 
decays Bg ^ — > 

On the other hand, along with the supersymmetric trilinear Yukawa interactions, 
the superpotential of the MSSM contains a bilinear SUSY term, ^HiH2, that couples the 
two Higgs-doublet superfields. This /i-term is phenomenologically necessary in order to 
break the chiral Peccei-Quinn (PQ) symmetry U{1)pq and to provide a mass for the as- 
yet experimentally unobserved higgsinos. The /i-parameter has dimensions of mass and 
is singlet under the Standard Model (SM) gauge group. It can therefore be driven by 
supergravity quantum effects to the mass scales Mqut or Mpianck- However, successful 
electroweak symmetry breaking in the MSSM implies that fi must be small, close to Msusy, 
namely at the scale at which SUSY is softly broken. This theoretical difficulty is often 
referred to as the //-problem of the MSSM [8,9]. 

The /i-problem may be resolved by allowing an effective fi parameter to be generated 
at the SUSY breaking scale, as the VEV of (the scalar component of) an additional gauge 
singlet Higgs superfield, S. The Higgs bilinear term of the MSSM superpotential is replaced 
with the trilinear coupling \SHiH2- In its simplest form, such a model re- introduces the 
undesirable ?7(1)pq symmetry of the superpotential, leading to experimentally excluded 
electroweak-scale axions. Several mechanisms for breaking the t/(l)pQ symmetry have 
been discussed in the literature [9] and lead to phenomenologically distinct models, such 
as the Next-to-Minimal Supersymmetric Standard Model (NMSSM) [10] and the Minimal 
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Non-minimal Supersymmetric Standard Model (MNSSM) [11]. 

An interesting limit for the NMSSM and the MNSSM is the so-called PQ symmetric 
limit. If f/(l)pQ is weakly broken, then the lightest CP-odd Higgs field Ai is predominantly 
singlet and very light, essentially becoming the pseudo-Goldstone boson of this symmetry. 
In this case, it is therefore natural to expect that a pseudoscalar particle Ai with a mass 
of a few GeV may have eluded detection so far. Production of such light, highly singlet 
pseudoscalars within the context of the NMSSM has been considered through the decays 
of SM-like Higgs fields [12], in associated production with charginos [13] and in rare decays 
of Upsilon mesons [14]. 

In the NMSSM and the MNSSM, the singlet Higgs fields have no tree-level couplings 
to the Standard Model (SM) fermions or gauge bosons. It has long been known that 1- 
loop threshold corrections can produce significant non-holomorphic Yukawa couplings in 
the MSSM at large tan/? [15-17]. Most recently, it has been realized [18] that analogous 
threshold corrections produce sizeable radiative Yukawa couplings even for the singlet Higgs 
bosons in minimal extensions of the MSSM. 

In this paper we consider the effects of FCNC Yukawa couplings of the singlet Higgs 
bosons within the MFV framework. Such effects have been considered in [19, 20] within 
the NMSSM in the limit of a light pseudoscalar singlet. The scope of our paper goes 
well beyond these studies, both analytically, as our effective Lagrangian approach includes 
a resummation of tan/5 enhanced terms, but also phenomenologically, as we discuss the 
implications of light singlets at or below the electroweak scale in both the NMSSM and the 
MNSSM. 

The structure of our paper is as follows: in Section 2 we present a manifestly gauge- 
and flavour- covariant effective Lagrangian framework for the calculation of FCNC vertices 
between the Higgs bosons and fermions. Section 3 summarises the relevant analytic results 
regarding FCNC 5-meson observables. Numerical results for various FCNC processes are 
presented in Section 4, for relevant scenarios in both the MNSSM and the NMSSM. Our 
conclusions are given in Section 5. 

2 Effective Lagrangian Formalism 

In this section we derive the form of the manifestly gauge- and fiavour-covariant effective 
Lagrangian for Higgs boson interactions with fermions. Our conventions and notations here 
closely follow those of [7]. The effective Lagrangian describing the down-type quark self- 
energy transition Q'j^ — > li^^ may be written in a gauge-symmetric and fiavour-covariant 
form as 

- /:fff[$i, $2, S] = d^^ (hM + Ah,[<|.i, $2, S])^, Q% + H.c. , (2.1) 
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where $1(2) are the scalar components of the Higgs doublet superfields giving masses to the 
down-type (up-type) quarks respectiveljE], S is the scalar component of the Higgs singlet 
superfield and is the 3x3 down-type Yukawa coupling matrix. In (12.11) the first term 
is the tree level contribution, whilst Ah^ is a 3 x 3 matrix which is a Coleman- Weinberg 
effective functional [21] of the background Higgs fields $1^2 and 5*. Observe that the 1- 
loop effective functional Ah^ has the same flavour- and gauge-transformation properties as 
hrf $|. Typical Feynman graphs that contribute to Ah^ to leading order in the Higgs fields 
$2 and S are displayed in Fig. [TJ In detail, the analytic form of Ah^ is 
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Here n = 4 — 2e is the number of dimensions in dimensional regularisation (DR), Iat is the 

X A^-dimensional identity matrix, r^'^'^ are the Pauli matrices, Pl(r) = ~ (+)75] 
the standard chirality projection operators and Cp is the quadratic Casimir invariant of 
QCD in the fundamental representation, i.e. Cp = 4/3. Moreover, we denote the U{1)y, 
SU{2)l and 5't/(3)c gauge couplings by (71,2,3 and the soft SUSY-breaking gluino mass 
by M3. Finally, and are 12 x 12- and 9 x 9-dimensional matrices which describe 
the squark and chargino-neutralino mass matrices, respectively, in the background of non- 
vanishing Higgs fields $1,2 and S. 

We express the 9x9 chargino-neutralino mass matrix Mc in the Weyl basis 

where Hu,d are SU{2)l doublets = (ht^Ki" 
Higgs field. In the NMSSM, Mc is given by 
/ 
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(2.3) 



"'^Here we adopt the convention for the Higgs doublets: = $2, Hd = «T2$|, where T2 is the usual 
Pauli matrix. 
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Figure 1: Gauge- and flavour-invariant self- energy graphs for the down-type quarks to 
leading order in the Higgs fields $2 o-nd S. 



where Mi 2 are the bino and wino soft SUSY-breaking masses and k is the trihnear singlet 
couphng of the NMSSM superpotential; the corresponding matrix for the MNSSM is ob- 
tained by setting k = 0. The 12 x 12 squark mass matrix is identical in both models. 
More explicitly, it is given by 



with 




In the above, h^j is the 3x3 up-type Yukawa coupling matrix, Mg are the 3x3 soft 
mass-squared matrices of the squarks and Sid^^ ~ ^d.u-^d,u 

are the corresponding 3x3 soft 
Yukawa matrices. Notice that (12.41) reduces to the MSSM result presented in [7], after 
identifying the effective /i-parameter as yU = — 4= Xvs, where vs = \pl (S) is the VEV of S. 
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The weak left- and right-handed quarks, denoted as m° ^, (f^ ^, are related to the 
respective mass eigenstates ul,r and dL,R by the unitary transformations 



R- 



(2.6) 



Here and U^'*^ are unitary 3x3 matrices, and V is the physical CKM-mixing matrix. 
The matrices and U^'*^ are determined by the mass renormalisation conditions 

;£ffj[$i,$2,S]) = -dnM^dL + H.c, (£:ff[$i,$2,^]) = -urMuUl + H.c, (2.7) 



where (...) indicates the vacuum expectation value of the enclosed expression. In (|27 
Mrf ji are the physical down- and up-quark mass matrices and C^f^ is the corresponding 
effective Lagrangian for up- type Yukawa sector [cf. fl2.ll) ]. Imposing the conditions 02.71) 
yields 



where 



V2 



V2 

V2 



(2.8) 



la + ^Uf (h;iAh„[$i,<|.2,5])uf 
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(2.9) 



Note that Ah„ is the respective 1-loop effective functional for the up-type quark Yukawa 
couplings. 

Considering now the general effective FCNC Lagrangian (12. ip for Higgs interactions 
to down-type quarks, we find that 
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where the component Higgs fields are given by 
$1,2 = 



1 



Ul 
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H.c, 
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(2.10) 



(2.11) 



In (12.1 op . A^^'^'^, A^^'^'"^ and A^^'^ are 3-by-3 matrices that are evaluated as 
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(2.12) 



where we have defined A^ = h^^ Ahd[^i, ^2, S] and suppressed the vanishing iso-doublet 
components on the LHS of (I2.12p . Since the singlet Higgs field does not contain a charged 
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component, the couplings of the charged Higgs boson are unchanged from those of the 
MSSM and we shall not consider them here any further. In the neutral sector, the physical 
Higgs boson mass eigenstates are related to the weak eigenstates by the orthogonal mixing 
matrices and O"^, such that 



/0i \ 
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(2.13) 



It is sometimes useful to parameterise by the CP-odd mixing angle 9a, where 



cos 6 a sin 6 a 
— sin 6 a cos 6 a 



(2.14) 



In terms of the mass eigenstates, the general effective Lagrangian (12.11) that includes 
the FCNC interactions of the neutral Higgs bosons with the down-type quarks reads: 
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where the Higgs couplings in the flavour basis = = = 13 are given by 
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(2.16) 
(2.17) 

(2.18) 
(2.19) 



In the above, we used the short-hand notation: tp = tan/3 = V2/V1, cp = cos/? and 
Sfs = sin p. 
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2.1 Single Higgs Insertion Approximation 



Our numerical results of Section H] include a full evaluation of and its derivatives using 
the matrix expression of (12.21) . We may obtain an understanding of our results by con- 
sidering an expansion of Ah^ in terms of $2 and S. We call the leading term of such an 
expansion, as represented by Fig. [H the single-Higgs-insertion (SHI) approximation. 

Within the SHI approximation, the tan /5-enhanced threshold corrections satisfy the 
simple relation 
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and /(x, z) is the 1-loop function, 

xy ln(x/?/) + yz ln(?// z) + xz \\i{z/x 



(2.20) 



(2.21) 



I{x,y,z) 



{x - y){y - z){x - z) 



(2.22) 



In writing the above, we have neglected the subdominant contributions coming from 
chargino- and neutralino-exchange graphs, which are weak-coupling constant suppressed. 

3 5-meson FCNC Observables 

In this section we review the analytic results relevant to FCNC 5-meson observables. Our 
conventions and discussion follow [6]. 

3.1 AMs, 

In the approximation of equal 5-meson lifetimes, the SM and SUSY contributions may be 
written separately as 
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(3.1) 



Although the SM predictions are consistent with the observed experimental values [22], the 
uncertainties are large and a non-negligible SUSY contribution is not excluded. 
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The SUSY contributions are given by 
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,(3.3) 



where DP stands for the Higgs-mediated double-penguin contributions. Here Fb^ is the 
weak decay constant and Bb^ is the so-called "Bag" parameter of the -6° meson respectively. 
The numerical factor r/^ is due to QCD corrections. We use the next-to-leading order QCD 
factors [5, 23] and hadronic matrix elements at the scale /i = 4.2 GeV: 
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The double-penguin Wilson coefficients in (13.21) are given by 
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where the relevant couplings 



^H,{Aj)bq 



are given in (12.161) . We neglect the 5-meson masses 



and Higgs boson widths in the denominators of the sums, except for the lightest Higgs mass 
eigenstates Hi and Ai, since we may have masses (9(1 GeV) in the case of a highly gauge 
singlet particle. In this limit we replace the effective propagators Mjj^f^j^^^ by their Breit- 
Wigner forms. Of the two relevant 1-loop contributions to (^B^\H^^'^\B^^ ^ the t — 
box contribution to 6*2^ may be given to a good approximation by [5] 
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We also include the chargino-stop box diagram which contributes to Cf^^, although its 
contribution remains sub dominant. 
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3.2 Bl^^fi+fi- 

Neglecting contributions proportional to the lighter quark masses nid^s, the effective Hamil- 
tonian for AB = 1 FCNC processes is given by 

H^/=' = -2V2GFVtbV,liCsOs + CpOp + CioOio) , (3.7) 

where 

_ _ 
Os = ~[^^b (qPRb) (jJ^fi) , 

_ _ 

Op = (QPRb) (/W75/i) , 

Oio = -^m, (qYPLb) (/i7M75/i) ■ (3.8) 
The Wilson coefficients Cs and Cp are given at large tan/? by 

and Cio = —4.221 denotes the leading SM contribution. Again we neglect the 5-meson 
masses and Higgs boson widths except for the contribution of the lightest Higgs scalar 
and pseudoscalar. The reduced scalar and pseudoscalar Higgs couplings to charged leptons 
^^HiA )Ti given in [18]. The leptonic vertex corrections can become important in the 
limit of a light, highly singlet Higgs boson, so we retain them in our numerical estimates. 

Using the effective /S.B = 1 Hamiltonian (13.71) . the branching ratio for i?^^ — > 
is given by [24] 
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(3.10) 

where q = d,s and r^^ is the total lifetime of the Bq meson. The form factors Fgp^ are 
given by 

Fl. = 4^1 F,.-!=L-C,P . n = -1f,.C,„ . (3.11) 
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4 Numerical Results 



In this section we provide numerical estimates of FCNC observables, within the contexts 
of both the MNSSM and NMSSM. We focus on those sectors of the models where the 
gauge-singlet Higgs fields are predicted to be light. We assume the framework of MFV so 
that all flavour changing effects are proportional to the CKM matrix V. This implies that 
the quark-sector soft SUSY-breaking terms are proportional to I3 [7], i.e. 

^l,u,D = Mlu,D I3, A„,, = Au,d I3 , (4.1) 

and similarly for the leptonic sector. We neglect the neutralino-mediated contributions to 
the FCNC Higgs boson couplings, which are found to be subdominant. 

Important constraints on FCNCs come from B Xgj experiments, which are sen- 
sitive to large stop-quark mixing due to a dominant chargino-stop quark loop. To avoid 
conflict with these constraints, we assume that Msusy is relatively heavy of order ~ 2 TeV. 
In calculating the Higgs couplings, we have used the following benchmark values for the 
soft SUSY-breaking parameters 

Ml = Ml = Ml = Ml = Ml = (1.7 TeV)2, 

Au = Ad = Ae = 2.0 TeV, (4.2) 
Ml = M2 = M3 = 2.0 TeV . 

In addition, we take /i = 140 GeV and t/^ = 50 throughout. The Higgs-mediated contri- 
bution to i? — > Xs'j proceeds through the charged Higgs boson and the predictions of the 
MSSM are not altered by the presence of the singlet Higgs bosons. We have made use of 
the public code CPsuperH [25] to check that all points considered here are consistent with 
the 2a experimental bounds on 5 — > Xg'j; the Higgs-mediated contribution was found to 
be subdominant in the region of interest. 

4.1 Electroweak-scale Higgs singlets in the MNSSM 

The renormalisable MNSSM superpotential is given by 

Wmnssm = Wvuk + XSH1H2 + tpS , (4.3) 

where VVVuk represents the Yukawa couplings of the MSSM superpotential. The tadpole 
parameter tp and its associated soft SUSY-breaking term ts are radiatively generated by 
supergravity quantum effects from Planck-suppressed non-renormalisable operators in the 
Kahler potential and superpotential. Discrete /^-symmetries, such as and Z^, that 
are imposed on the theory postpones the appearance of these operators beyond the 5-loop 
level, such that they are naturally suppressed of the order of Msusy within a perturbative 



11 



framework of supergravity (SUGRA). For a detailed proof of the argument, see the first 
reference of [11]. 

The tree-level Higgs sector of the MNSSM may be described by the six parameters 

tp, Ma, /i, A, m?2, — , (4.4) 

where Ma is the would-be MSSM pseudoscalar Higgs mass and = ^tp- In our numerical 
results we also include the dominant radiative corrections due to both (s)top and (s)bottom 
loops, which are necessary to raise the mass of the SM-like Higgs boson above the LEP II 
constraint of 114 GeV. 

Our interest here is in the effects of light gauge-singlet Higgs bosons, with masses 
at or below the electroweak scale. Such a scenario may be realised by taking Ma to be 
large, so that the heavy Higgs doublet fields effectively decouple, and by choosing such 
that the mixing between the light CP-even Higgs bosons is suppressed. In the numerical 
estimates presented here, we take 

M„ = 1.5TeV, = (1.0 TeV)^ /i = 140 GeV, tan/3 = 50, A = 0.3 . (4.5) 

This corresponds to a Higgs pseudoscalar mixing angle of cos 6*^ ~ 0.17. The scalar Higgs 
singlet mixes almost exclusively with the 0i doublet, the non-singlet fraction is approx- 
imately equal for both the lightest Higgs states Hi and Ai. These light particles are 
not constrained by the direct search limits from LEP II and the Tevatron due to their 
highly gauge-singlet nature, with gHi,2ZZi QhiAiZ < 0.05. The mass scale of the singlet 
Higgs bosons are approximately given by ^Jxtsjjl. Notice that, in the absence of strong 
singlet-doublet mixing effects, the singlet Higgs scalar and pseudoscalar are approximately 
degenerate due to the tree- level mass-sum rule [11]: 

+ + = Ml + + m% . (4.6) 

Including dominant radiative corrections we find the masses of the remaining Higgs bosons 
in this scenario to be 

= 147 GeV, m^g = 1.52 TeV, m^^ = 1.52 TeV, (4.7) 

so that the H2 Higgs boson, which has SM-like couplings to the EW gauge bosons, is well 
above the direct search limit set by LEP II. 

4.1.1 Effects on AMb, 

In Fig. [2] we show the SUSY contribution to AMb^ in units of pS , cLS cL function of the 
lightest pseudoscalar mass. The upper curve fully includes the radiative Yukawa couplings 
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Figure 2: The SUSY contribution to AMb^ in units ofps~^ as a function of the mass of the 
lightest pseudoscalar niA^ in the MNSSM. All parameters are taken as in Iji4-^ ^'^^ 
The horizontal lines show the currently measured value along with the SM and MSSM 
predictions for corresponding values of Ma and tan/5. 



of the Higgs singlet fields whilst the lower curve neglects these corrections. We observe 
that the SUSY contribution can be larger than the currently observed value at large tan (3 
if the singlet Higgs bosons are light. For moderate values of m^i, the effects of the direct 
singlet Yukawa coupling remain significant, enhancing the value of AMb^ by around 25%. 

The dominant contribution to AMb^ for light singlet Higgs bosons is due to the Wilson 
coefficient Considering the SHI approximation for the reduced couplings SHi(Ai)bq' 

we see that the off-diagonal terms are given by 

g^,5, = (t,Of, - ^Oa^l) VtR-^V , (4.8) 
= ^ (t.Ofi + ^O^l) VtR-^V , (4.9) 
where we have neglected the small contributions of the $2 component next to those of $1. 
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Figure 3: The dominant Wilson coefficients contributing to AMb^ as a function of the 
mass of the lightest pseudoscalar in the MNSSM. The dashed curve shows the Wilson 
coefficient neglecting the threshold corrections for the gauge singlet Higgs boson. 



We observe that the FCNC couphngs of the hghtest scalar and pseudoscalar are of equal 
magnitude, if threshold corrections to the singlet-Higgs Yukawa couplings are neglected 
and Cf( ~ In this case, there is a cancellation between the Hi- and the Ai-mediated 
contributions to An analogous cancellation between the heavy Higgs bosons H 

and A is known to take place in the MSSM as well, as a consequence of the PQ symmetry 
which forbids at the tree-level the dominant AB = 2 operators, [bRdiy and (bRSL)'^, 
associated with (7^^^(°p) (^gee, e.g. [5]). 

For very light singlets this cancellation is dominantly broken by the mass splitting 
between Hi and Ai, typically of the order rriHi ~ ttlji-^ + 2 GeV. For larger masses, this 
splitting is negligible and the dominant breaking is instead due to threshold effects on the 
singlet-Higgs Yukawa couplings, which are found to contribute to the reduced couplings 
with opposite sign, due to the mixing matrices. The dominant Wilson coefficients for 
AMb^ are plotted in Fig. [3l We do not plot (^^^^(^p)^ which exhibits a behaviour similar 
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Figure 4: The SUSY contribution to ^Mb^ in units of ps^^ as a function of the mass of 
the lightest pseudoscalar rriAi in the MNSSM. The line conventions are as for Fig. [H All 
parameters are taken as in {j-'^ o.i^'d ( f^.^l j. The upper horizontal line shows the currently 
measured value and the lower horizontal line shows the MSSM prediction for corresponding 
values of Ma and tan/5. We do not show the SM prediction as the central value is close to 
the experimentally observed splitting. 



to (j^^^^^^^ ^ which is suppressed by a factor of {ms/mb)'^. 

In Fig. Hlwe show the corresponding plot for the contribution to AMb^- Again we 
observe that the SUSY contribution can exceed the currently observed value, with the 
limit on m^^ at the same level as that due to AMb^- The effects of the singlet Yukawa 
coupling are more pronounced at larger values of m^^ than in the case of AMb^, as the 
enhanced Wilson coefficient remains dominant over the two-Higgs-doublet model 

contribution (^^r(2Hdm) order of magnitude here even for niAi ~ GeV. 
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Figure 5: The branching ratio B{Bs — * as a function of the mass of the lightest 

pseudoscalar niA^ in the MNSSM. All parameters are taken as in Iji4-^ ^'^d i4-^ - The 
horizontal lines show the SM [26] and MSSM predictions and the current experimental 
upper limit at 90% C.L. 



4.1.2 Effects on B^^ 

Figure [5] shows the branching ratio B{Bs function of m^^ . As expected, the 

prediction varies strongly with the mass of the hght gauge singlet Higgs bosons due to the 
^Ht(Ai) dependence of the form factors -F|(p)- The prediction exceeds the current bounds 
for Higgs singlet masses below around 50 GeV. Since the branching ratio in (13.101) depends 
only on the absolute values of the form factors, there is no cancellation between the scalar 
and pseudoscalar contributions. At low values of m^^, including the threshold corrections 
leads to a slight increase in Cp, which is compensated for by a slight decrease in Cs and the 
overall effect is negligible. At higher masses, m^^^60 GeV, there is cancellation between the 
contributions due to Cp and Ciq. In this region the singlet threshold corrections produce 
a noticeable shift. 
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4.2 Light singlet Higgs pseudoscalar in the NMSSM 



The NMSSM superpotential is given by 



Wnmssm = VVvuk + XSH1H2 + —S^ . 



(4.10) 



The corresponding soft SUSY-breaking terms are given by 



In the hmit A 




symmetry beyond t/(l)pQ. This symmetry results from an /^-symmetry, denoted by [/(1)r, 
which is exact in the absence of soft SUSY-breaking trihnear couphngs. Recall that the 
superpotential of the NMSSM is charged under U(1)r and so it differs from f/(l)pQ. This 
U{l)n is spontaneously broken when 5" acquires a VEV, so that for small values of Ax^k the 
singlet Higgs pseudoscalar is again a pseudo-Goldstone boson and is naturally expected to 
be much lighter than the other Higgs fields [11,27]. 

In order to examine FCNC observables within such a scenario, we take the CP-odd 
mixing angle 6a and the mass of the lightest Higgs pseudoscalar m^^ to be free parameters 
in place of Ax ,^, along with the couplings A and k. We use the following benchmark values 
throughout this section; 



which are found to lead to the following masses for the remaining Higgs bosons, 

ruH, = 135 GeV, tuh^ = 357 GeV, mn., = 1.14 TeV, tua^ = 1.14 TeV. (4.13) 
The lightest CP-even Higgs in this scenario has SM-like couplings to the EW gauge bosons. 

For general values of n, the CP-odd and CP-even singlets of the NMSSM are not 
constrained to have degenerate masses. In the scenario considered here, there is only one 
light Higgs particle, the pseudo-Goldstone boson Ai. Due to this, there is no cancellation 
between the dominant Higgs field contributions to the Wilson coefficients as described 
in Section 14.11 for the MNSSM. This forces us to take a smaller singlet-doublet mixing 
parameter cos^^ in order to find phenomenologically acceptable results at large values of 
tan p. Note that the value of cos 6'^ used here is only ~ 10% of that found for the MNSSM 
scenario considered previously and that the contributions to g^^^;^ from singlet-doublet 
mixing and direct threshold corrections are therefore comparable in magnitude. 

4.2.1 Effects on AMb^ 

In Fig. E] we show the SUSY contribution to AMb^ as a function of the lightest Higgs 
pseudoscalar mass. In this scenario the SUSY contribution also exceeds the currently 
measured value of AMb^ for the lightest values of rriAi at large values of tan p. 



A = 0.4, K = -0.5, coseA = 0.018, /i= 140 GeV, tan/3 = 50 , 



(4.12) 
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Figure 6: The SUSY contribution to AMb^ in units of ps^^ as a function of the mass 
of the lightest pseudoscalar niA^ in the NMSSM. The line conventions are as for Fig. 
All parameters are taken as in { 4-2 } and The horizontal lines show the currently 

measured value along with the SM and MSSM predictions for corresponding values of Ma 
and tan/3. 



We remark that the contribution to AMb^ in this scenario of the NMSSM is rather 
smaller than in the MNSSM scenario discussed in Section HTTl The leading contribution in 
both cases is due to the Wilson coefficient \^ discussed above, in the NMSSM 

there is no cancellation between dominant CP-odd and CP-even Higgs fields and so we 
have taken the FCNC Yukawa couplings to be an order of magnitude smaller than those 
considered in the previous section. As a result of this, the contribution of the singlet 
pseudoscalar become negligible here for mAi~25 GeV and we recover the MSSM prediction. 

Figure [7] shows the corresponding predictions for AMb^- As in the MNSSM scenario 
considered in Section 14.11 the bounds from — mixing are at the same level as those 
from the Bg mesons, although the relative importance of the singlet threshold corrections 
at larger masses is greater due to the continuing dominance of (j^^^^^'> over C^^^^^\ In 
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Figure 7: The SUSY contribution to ^Mb^ in units of ps^^ as a function of the mass of 
the lightest pseudoscalar rriAi in the NMSSM. The line conventions are as for Fig. [H All 
parameters are taken as in { 4-2 } and ^^A^. The upper horizontal line shows the currently 
measured value and the lower horizontal line shows the MSSM prediction for corresponding 
values of Ma and tan/5. We do not show the SM prediction as the central value is close to 
the experimentally observed splitting. 



particular, the contributions of the singlet pseudoscalar are not neghgible at the same low 
masses as for AMb^- 

4.2.2 Effects on 5° 

Figure E] shows the branching ratio B{Bs as a function of the lightest pseudoscalar 

mass in the NMSSM. For small values of m^^ the prediction exceeds the current 90% 
confidence limits at large tan/3. For a heavier singlet Higgs pseudoscalar, we observe 
that inclusion of the singlet threshold corrections leads to a significant reduction in the 
branching ratio, which becomes suppressed by more than an order of magnitude compared 
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Figure 8: The branching ratio B{Bs cls a function of the mass of the lightest 

pseudoscalar, m^i in the NMSSM. The solid line includes threshold corrections for the 
gauge singlet Higgs bosons, the dashed curve neglects these. All parameters are taken as 
in lli4-2^ and l[4-12 ). The dot-dashed curve shows the NMSSM prediction, including all 
threshold corrections, but with = —2 TeV. The horizontal lines show the SM [26] and 
the MSSM predictions and the current experimental upper limit at 90% C.L. respectively. 



to the SM prediction. In this mass range the contribution of the hght pseudoscalar interferes 
destructively with the SM-like diagrams. It should be noted that such a large suppression is 
not possible in either the MSSM or MNSSM, since in both these models tree-level mass-sum 
rules prevent the appearance of isolated CP-odd Higgs bosons. 

The leading effects of the neutral Higgs bosons on B{Bs f^^fJ'^) in both the MSSM 
and MNSSM may be understood by including only the contributions due to the dominant 
scalar /pseudoscalar pair (i.e. H,A in the MSSM or Hi,Ai in the above considered limit 
of the MNSSM) and by assuming equal masses and effective Yukawa couplings (this is 
accurate to the order of a few percent). In this approximation, it can be shown that the 
absolute minimum of B{Bs /i^A*") in both these models is 1/2 x B{Bs /i"*'/^^)sM- 
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When such a cancellation is effective in the NMSSM, the absolute minimum of 
B{Bs yU^/i") is due to the contribution of the heavy CP-even Higgs doublet, H3 ~ (pi. 
The exact value of B{Bs at the minimum scales as Mjj'^ and the branching 

ratio will be further suppressed as and A2 become heavier. As can be seen from the 
dot-dashed curve in Fig. [H the suppression is found to vanish when ^A^ flips sign, i.e. for 
Au = —2 TeV, since the leading SM and SUSY contributions interfere constructively in 
this case. 

5 Conclusions 

We have shown that singlet Higgs bosons generically present in extensions of the MSSM can 
contribute significantly to FCNC observables at large tan j3, both through their tree-level 
mixing with the Higgs doublet fields and through 1-loop quantum effects. For very light 
Higgs singlets, the contributions may exceed even by one order of magnitude the current 
experimental limits on FCNC processes. This allows us to place severe constraints on the 
mass of the lightest scalar state in the theory at large values of tan j3. 

In the MNSSM, the scalar and pseudoscalar singlet Higgs fields are constrained to 
be approximately degenerate by a tree-level mass sum rule [11]. Their contributions to B- 
meson mixing through the dominant Wilson coefficient Cf^^ interfere destructively, analo- 
gously to the contributions of the MSSM heavy Higgs bosons. At moderate to large values 
of the common Higgs singlet mass, the direct coupling plays an important role in lifting 
this cancellation, leading to a large enhancement of Cf^^. Although cancellations between 
the pseudoscalar and SM-like contributions to the branching ratio Bg f^^f^~ are possible, 
the quasi-degenerate scalar Higgs singlet produces a screening effect. Suppression of the 
branching ratio below the SM prediction is limited to at most a factor of ~ 1/2 at leading 
order. 

The situation in the NMSSM is different. The CP-even singlet remains heavy close 
to the i?-symmetric limit. The contribution of the singlet pseudoscalar can therefore easily 
exceed the experimental limits at large values of tan (3, unless the double-singlet mixing is 
constrained to be small. Due to the absence of an accompanying CP-even Higgs boson, the 
cancellation between SUSY and SM-like contribution to the branching ratio B{Bs f^^f^~) 
can be highly efficient in this model, where suppression can be more than one order of 
magnitude relative to the SM prediction. Such an effect is not possible within the MSSM 
at large tan (3. Possible discovery of SUSY partners at the LHC along with non-observation 
of the decay B^ l^'^ would strongly point towards the NMSSM and searches for light 
pseudoscalar particles should become a high priority, should such a situation present itself 
in the near future. 
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The present work was carried out within the framework of minimal flavour violation. 
However, within a general renormalisable framework of the models we have been studying 
here, new sources of CP violation could be present amongst the soft SUSY-breaking terms. 
This is expected to lead to new CP-violating threshold corrections for both the singlet and 
doublet Higgs fields. Additionally, both the NMSSM and the MNSSM exhibit explicit CP 
violation in the tree-level Higgs potential. These effects may lead to significant enhance- 
ments of the signals considered here, along with contributions to other observables, such 
as electric dipole moments. We plan to report progress on these issues in the near future. 
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A The Chargino-Neutralino Propagator Matrix 



In ( I2.2p . we need to evaluate the different left- and right-handed chiral components of the 
chargino-neutralino propagator matrix: 

To do so, it proves more convenient to use the chiral representation for flA.ip and write 
Ac(^) as 

Ac(« = f-^e^^ Y . (A.2) 

In the above, we have defined k = k^a^ and k = /c^a^, where = (I2, cr), = (I2, — cr) 
and (T are the usual Pauli matrices. Note that kk = kk = k^l2- 

Our calculational task is equivalent to finding the inverse of an invertible, block 
matrix M, 

where A, B, C and D are m x m-, m x n-, n x m- and x n-dimensional matrices, 
respectively. To this end, we denote the inverse of M by 

Correspondingly, a, /3, 7 and 5 are matrices of dimensionality mxm, mxn, nxm and 
n X n. We may now use the condition MM~^ = to write the four constraints, 

Aa + B-f = 1^ , (A.5) 

Ca + D^ = 0„xn^, (A.6) 

AP + B6 = 0„,xn, (A.7) 

C13 + D5 = In . (A.8) 

Conditions (1A.6P and (1A.7P imply that 

7 = -D-^Ca , (A.9) 
f3 = -A^^B6 . (A.IO) 

Combining these with (1A.5I1 and (1A.8I) . we may derive 

a = (^A-BD^^Cy^ , (A.ll) 
5 = [p-CA-^B)'^ , (A. 12) 
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so that the inverse matrix M ^ may be written as 

\ -D~^C{A- BD-^Cy^ {D-CA~^By^ j ' ^ ' ' 

One should observe here that the expression is not unique. Specifically, we could have 
used M~^M = Im+n to obtain equations analogous to flA.5l) - flA.8l) and express in a 
form which, although less compact, relies on the inverses of two combinations of submatrices 
only, e.g. D-^ and {A - BD-^Cy . Such expressions are more efficiently implemented in 
numerical routines. However, all these alternative forms are equivalent to each other. 

We now return to ^c{¥) i'^ (lA.2p . Using (lA.lSp . the chargino-neutralino propagator 
matrix may now be written down as 



^ M[, (Pig - McMU ' I2 (Pig - mJ,Mc) ' k 



y yPlg-McM^) 'k Mc(Plg-M^Mc) 'I2 



(A. 14) 



where we have used the fact that the only non-zero commutator amongst the submatrices is 
[A, D] which simplifies the above expression. From (]A.14p . we may now obtain the different 
block elements, by appropriately acting with left- and right-handed projection operators 
Pl.r- Foi' example, for the expression that occurs in (12.21) . we find that 

PL^CmPL = MUPlg-McMU"'l2 , (A.15) 



where we have suppressed the vanishing elements on the RHS of ( ]A.15|) . 



Using flA.lSp . we may numerically integrate expressions such as (12. 2p without any 
need to first diagonalise the mass matrices appearing in the propagators. At each sample 
point in the integral, the propagator matrices are numerically inverted and the matrix 
element relevant to the process is selected, i.e. the propagator matrix element representing 
the transition between the required electroweak eigenstates. In this way, we avoid the 
introduction of cumbersome mixing matrices for the squarks and charginos in the analytic 
expressions of Section O 

Furthermore, by working with all propagator matrices in the electroweak basis, the 
calculation of the derivatives appearing in (I2.12p is greatly simplified. The expression for 
each derivate of the propagators in the matrix form contains no more than two terms. By 
contrast, if we were to diagonalise the mass matrices, we would find expressions containing 
three terms for each mass eigenstate (coming from applying the Leibniz rule to the product 
of Higgs-dependant masses with initial and final mixing matrices), many of which typically 
require further expansion before one could cast them into a useful compact form. 
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